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It is shown that quantitative information on spatial correlations in a system of polarizable particles can be 
extracted directly from its experimentally measurable optical spectra. For a collection of metallic nanoparticles 
(NPs), it is demonstrated that the degree of short-range correlation in NP's positions can be evaluated by an 
appropriate numerical analysis of the extinction spectrum in the surface plasmon resonance region, given the 
polarizability of an individual NR The spectrum analysis consists in the evaluation of a single number, which 
is the derivative of the ensemble response function in the vanishing polarizability limit, to be compared to 
pre-calculated values for a model NP system with given density and correlation parameters. 

PACS numbers: 78.67.-n, 81.07.-b 



Introduction. The optical properties of assemblies of al- 
most identical polarizable particles, randomly distributed in 
space and coupled by electromagnetic interactions have at- 
tracted attention for a long time 1 1 , 2], since they are of inter- 
est for fields ranging from interstellar dust [ 3 ] to electrorhe- 
ological fluids [4 ] to biosensors and single molecule detec- 
tors (5] |6) • Metallic nanoparticles (NPs) dispersed in a di- 
electric matrix are a practically important example of such as- 
semblies (7). These nanocomposite materials have been pro- 
posed for environment sensing 00 [9), tailoring color of func- 
tional coatings [10], enhancement of the non-linear response 
(IUE1, energy harvesting in solar cells fT3l and Raman scat- 
tering by molecules fl4l . amongst other applications of the 
nanophotonics. Theoretical description of the optical proper- 
ties of a composite containing NPs becomes particularly inter- 
esting (and complex) when the polarizable particles are non- 
randomly distributed in space fT5llT7lL e.g. aggregate into 
clusters, possibly of a fractal dimension |[l6l[T8]|. Note that 
NP aggregation is precisely one of the mechanisms proposed 
for sensing (6j |9) ; also it can enhance the detection capability 
of surface-enhanced spectroscopy methods fT9l . 

Is it possible to extract information concerning the spa- 
tial correlations in the particles' positions from the exper- 
imentally measurable optical response of the system with- 
out performing tedious fitting of the spectra using a phys- 
ically founded numerical model? In the present work, we 
address this fundamental question by studying ensembles of 
nanospheres described by so called coupled dipole equations 
(CDE) (21Q21IT71. CDE are linear algebraic equations relat- 
ing the polarization of each particle to the external electric 
field. The dipole-dipole interactions between the particles are 
represented by a second rank tensor and CDE contain a spec- 
tral variable and some parameters depending on the material 
properties. The purpose of this work is to study the relation 
between the spectral-dependent solutions of the CDE and the 
type and parameters of the particle-particle correlation func- 
tion in real space. It may seem desperate trying to solve the 
inverse problem in order to obtain the particle-particle corre- 



lation function from an experimentally measurable spectrum 
of some quantity, such as absorption or extinction, because 
the spatial correlations affect the collective optical properties 
in a rather indirect way. However, one may hope to solve a 
more modest but still relevant problem consisting in the use of 
the optical spectroscopy data in order to distinguish between 
some classes of spatial distributions of the particles, e.g. be- 
tween a completely random distribution and a short-range cor- 
related one where the particles form closely- spaced pairs. Our 
calculations are performed for a collection of metallic NPs 
(with the dielectric function described by the Drude model), 
distributed over a square or a cubic lattice. We will show that 
it is possible to propose a numerical criterium that allows for 
conclusions concerning the degree of short-range correlation 
in the particles' positions using the (experimentally measur- 
able) extinction spectrum of the NP system. The central result 
consists in the calibration of this criterium as a function of 
two arguments, NP concentration and short-range correlation 
parameter, thus making it possible to classify an ensemble of 
NPs according to their collective optical response. Therefore 
one can hope that, at least, some classes of non-random distri- 
butions of the NPs can be distinguished by means of a rather 
simple analysis of their extinction spectra. 

Coupled Dipole Equations. First, we present the formalism 
used to describe the electromagnetic (EM) interactions be- 
tween polarizable particles that randomly occupy some frac- 
tion of the sites of a regular lattice. If the characteristic wave- 
lengths of the EM interactions involved are much larger than 
the size of each NP and the distance between them, the CDE 
system can be written as 12: 

di = a ^E () + ^T /y d y j , (1) 

where i and j run over all N lattice sites, is the dipole 
moment located at site i, Eo is an external electric field 
that will be taken as uniform and directed along the z axis 
in 3D, E = (0, 0, Eq), where E is a constant. In 
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FIG. 1 : (Colour online) Ensemble response function spectra, normal- 
ized to the maximum value, r = c/)(ao(uj))/ max[0]), for random 
(£ = 0, full curve) and strongly clustered (£ = 20, dashed curve) 
distributions of Au NPs on square lattice (c = 0.2, a /I — 0.25). 
SPR mode of a single NP is located at 2.37 eV, the peak at 2.31 eV is 
due to pairs of NPs occupying sites that are nearest neighbours along 
the direction of the applied field. The inset shows the hodograph of 
the function ao for a Drude particle with oj p — 6.91 eV, Coo = 1.53 
(bulk Au parameters), T p calculated as in Ref. 1 10 1 and en — 5.75. 



Eq. |IJ) Tij is the dipole-dipole interaction tensor, = 
R^ 3 (3n a np — 6 a p), where Rij is the distance between lat- 
tice sites i and j, a and f3 denote Cartesian components and 
n is a unitary vector in the direction from site i to site j. 
Let us express it in units of the lattice constant, I, so 
is dimensionless and corresponds to a quasistatic approxi- 
mation, in which is a real and symmetric matrix. The 
NP's polarizability, ao, is a complex function of frequency, 
uj. It is an experimentally controllable parameter. If the 
dipoles are spherical metallic (nano) particles, then we have 
ao(cj) = a 3 (e(uj) — 1)/(c(uj) + 2), and, within the classical 
Drude model [7], e(uj) = (e^/e^) (l - uj%/(uj(uj + iT p ))). 
Here, a is the NP radius, eh is the dielectric constant of the 
host matrix, and e^, uj p and T p are the parameters of the 
plasma resonance in the NP material. 

Let V be the dimension of the space (V = 3 if we are deal- 
ing with a 3D problem and V = 2 in the 2D case). Denote by 
\p) a vector with VcN components (c is the fraction of lattice 
sites occuppied by the particles). Let / be the set of indices 
corresponding to the occupied sites. Then we have \p) = 

where the index is included if and only if it belongs to /. 
We define the dimensionless polarization vector components, 
Pi a , aspia = d ia /(Eoao), i e I. The constant vector \e) has 
components given by = e a = (0, 1), i G I,ifV = 2, 
and, ei a = e a = (0, 0, 1), i G /, if V = 3. By M we denote 
the VcN x DcTV-matrix obtained from the VN x VN matrix 
formed by blocks T^, i, j = 1, N, eliminating all zero rows 
and columns. This is a dimensionless symmetric matrix with 
zero diagonal elements. 



Now we can formulate the problem in the matrix form. 
Equation ([T]) is rewritten as 



\p) = \e) + a M\p) . 



(2) 



Its solution has the form \p) = (J — aoM) _1 |e) = 
G\e), where G is Green's function. If we denote the 
eigenvalues of M by (k = 1,2..., VcN) and the 
corresponding eigenvectors by \uk), then G (ao(u)) = 
Ylkli \ u k)(uk\(l — c^o^/c) -1 - It is a VcN x VcN matrix 
depending on the spectral variable. Although normally the 
EM frequency (uj) is used as such, we may also consider the 
polarizability (ao) as a complex spectral variable. 

Most observable quantities of interest can be expressed via 
the Green function. One of these, the extinction cross section 
characterizing the processes of both absorption and scattering 
of light by the particles, is given by l20l : 
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where A denotes the wavelength, F(S n ) is the distribution 
function of N± possible configurations of the A^-dimensional 
vector, 5 =(..., 5i,. . .), that satisfy the condition, \S\ = cN. 
The number N\ is given by Ni = N\/((cN)\[(l - c)N]\), 

T,n=i F ($n) = 1, and Yln=i S i F ($n) = c, for any site 
i. The scalar product(e|p) n can be expressed through the 
Green function of the corresponding configuration of NPs: 
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e ) . Thus, the calculation of the optical 



response of a system of coupled dipoles with known F(S n ) 
is a straightforward (although possibly computationally hard) 
task. Two examples of calculated response spectra are shown 
in Fig. [T] NP pairs oriented parallel to the applied field pro- 
duce the low energy peak in the strong clustering case, as it 
has been observed experimentally I2H . 

Ensemble Response Function. Our purpose is to define 
a quantity that is connected with the distribution function, 
F(S n ), and, at the same time, directly related to the spectral 
properties of G. Let us consider the function 
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The values of (j)(a), which we shall call ensemble response 
function (ERF), are known for a = ao(u) [i.e. on the hodo- 
graph shown in the inset of Fig. [TJ, since they can be obtained 
from experimental data for the extinction cross-section. 

Let uj* be such that Re(ao)(^*) = 0. Set a* = ao(uj*). 
The function (/>(•) is analytic in a disk {a | \a — a*\ < 
\/a% + A 2 ; A>0, A<<1} containing zero. In a small 
neighbourhood of zero we have: 

0(a) = £ [<e \e) F(5 n ) + a (e \m M | e) F(5 n ) + . . . 

n=l 

= Nc + aY, Kij e « e ^f + • • -(5) 
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We are interested in cases where the (random) occupation 
of the lattice sites is correlated. This is described by the 
particle-particle correlation matrix defined as follows: 
Kij = Yln=i^jF(5n)' Then it follows from Eq. ^ that 
the ERF derivative taken at a = is directly related to the 
correlation matrix, 

4>'(0) = J2K ij J2e<*e l 3T?f . (6) 

i,j aj3 

The ERF derivative value ([6]), with an appropriate normal- 
ization, is the proposed criterium for the evaluation of the de- 
gree of particle-particle correlations in the NP ensemble. Ob- 
viously, it is not possible to extract the correlation matrix el- 
ements from the optical response spectrum, not even the pair 
correlation function, k(r) = ^Z^ =1 J2(i>j- R tj =r) K ij ( wnere 
r is the distance between the 1-st, 2-nd, etc nearest neighbors 
in the lattice). This single number, obtainable from an exper- 
imentally measured spectrum, allows one to classify the NP 
ensemble as e.g. a strongly aggregated or a gas-like uncorre- 
cted system. Of course, one needs calibration of the possi- 
ble values of this quantity, which can be obtained by solving 
CDEs for a set of model systems where the particles are dis- 
tributed according to predefined rules. Then the ERF deriva- 
tive value obtained from the spectrum of the real system will 
show its similarity to one of these models. This is much less 
time consuming than the fitting of the experimental spectra. 

So, technically the problem is reduced to the evaluation of 
0'(O) from the known values of ^{olq{uj)). It can be done by 
means of an appropriate extrapolation procedure that will be 
used for both real and model spectra. For real spectra, this is 
a hard ill-posed problem (note that, according to Eq. ([3]), the 
observable quantity that will allow us to compare simulation 
results with reality is proportional to 9ra(ao0(ao))). How- 
ever, since the function is a rational function with poles of 
(j)(a) concentrated in the real axis, the Tikhonov regulariza- 
tion method (22l|23] gives good results. 

Examples. For a 3D lattice and z axis chosen along the 
direction of the external field, 0'(O) = K ij T if- If the 
lattice is cubic and the particle distribution is isotropic, we 
can write in the continuum approximation: 

0'(O) = J k(r)^^^-r 2 dQdr . 

The integration yields zero whatever the pair correlation func- 
tion is, since it depends only on the modulus of the radius- 
vector. Therefore, for an isotropic 3D ensemble of spherical 
particles, the ERF derivative is equal to zero. This is an im- 
portant result, although it excludes the possibility of probing 
the correlations in this case. 

Let us consider a square lattice, which can be seen as an 
ultimate anisotropic limit of the 3D case. Nearly 2D struc- 
tures can be built from metal NPs grown by colloidal chem- 
istry techniques l5l [T9l [24^261 . Aggregation of NPs in such 
layers can occur as a result of some chemical treatment l26l . 
We shall limit ourselves by the case of short-range correla- 
tions that can be described by a single short-range order (SRO) 



parameter, £. If one particle is located at a lattice site ii, the 
relative occupation probabilities by another particle for the re- 
maining N — 1 sites (i ^ i\) are given by: Pi/i x = c + A^£, 
where A iil = 1 if R iil = 1, and A^ = otherwise. If 
£ > 0, this distribution corresponds to NP clustering. Using 
the Monte Carlo (MC) method, we generated a sufficiently 
large number of samples containing cN particles distributed 
over N = N 2 sites of a square lattice, and calculated the ERF 
spectra [by solving Eqs. |T])] and the correlation matrix (Kij), 
for given values of £ and c. Two examples of the calculated 
spectra are shown in Fig. [T] Since the eigenmodes of the 
system (for each MC sample) are known in this case, the nor- 
malized ERF derivative, t'(0), could be calculated explicitly 
and MC-averaged, as a function of c and £. As we established 
empirically, the values of t'(0) become almost independent 
of the lattice size for N x > 20. Of course, the results depend 
on the single particle polarizability and its definition requires 
the knowledge of several essential parameters to be obtained 
from the experiment. 

The calculated results are shown in Fig. [2] First, in order 
to check the validity of the extrapolation procedure described 
above, it was applied to the calculated ERF spectra treated as 
"experimental" ones in this case. We verified that the SRO 
parameter obtained by applying this procedure coincides with 
the true one (used in the MC method), with a relative preci- 
sion better than 1%. We also applied the extrapolation proce- 
dure with regularization to real experimental spectra of the ex- 
tinction coefficient [ 27 ] measured for thin colloidal gold films 
|25l . Although their colloidal layers were not exactly 2D NP 
systems (they consisted of approximately two NP monolay- 
ers), this should not be a serious problem, because we use 
normalized spectra. The mapping of this system to our lattice 
model is as follows, (i) Assuming that NPs can touch, we take 
I = 2a atp, where a atp is the "physical NP radius" determined 
by means of transmission electron microscopy [25 ]. Note that 
a atp > a since the NP's Au core was covered with a protect- 
ing shell, (ii) The Drude model parameters were either taken 
as for bulk gold (oj p , e^) or evaluated using the experimen- 
tal spectra of Ref. [25 ]. Namely, from the broadening of the 
extinction spectrum of NPs in colloidal solution we obtained 
T p « 0.45 eV and from the peak position in the spectrum of 
the most diluted film we found the "host dielectric constant", 
€h ~ 5.8. Note that the latter is an effective parameter, since 
the particles located at the substrate/air interface and in each 
NP its metallic core is surrounded by its shell. Using a known 
semi-empirical relation [24 ] between T p and the metallic par- 
ticle radius, we can estimate the core radius as a « 3 nm. 
Therefore, we used a/ 1 = 0.25 when computing the calibra- 
tion data presented in Fig. [2] 

The values of the "clustering" criterium for three NP films 
studied in Ref. 25 were calculated by roughly approximating 
the experimental spectra with Lorentzians and using the ex- 
trapolation procedure with regularization. They are shown in 
Fig. [2] defining "compatibility curves" (c, £). If c is known 
experimentally, the clustering parameter is determined. For 
the present example, the fraction of occupied sites can be 
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FIG. 2: (Colour online) Normalized ERF derivative at a — plotted 
versus NP concentration and SRO parameter for a collection of gold 
NPs on square lattice. Level curves for three experimental values of 
c are shown. The parameters are the same as in Fig. [T] 

related to the NP volume fraction / (provided in Ref. l25k 
by c = 6//7T. This way we obtained the values of c = 
0.11; 0.18; 0.24 for the three films, respectively. Thus, we 
can conclude from Fig. [2] that there was very little cluster- 
ing in these films, excepting perhaps the lowest NP concen- 
tration one (£ = 0.73), as qualitatively confirmed by their 
microscopy images (25). 

In summary, we described a mathematical procedure allow- 
ing for the evaluation of the correlation parameter, which is a 
characteristic of short-range pairing of the polarizable par- 
ticles, from the spectrum of an experimentally measurable 
quantity 3ra(ao<K<2o)) (normalized to its maximum value). 
To demonstrate the correctness of this procedure, we numeri- 
cally generated "experimental" spectra of a model 2D system 
of Drude particles assuming the nearest-neighbour pairing 
with a predefined SRO parameter (£) and obtained its value 
with high precision using the extrapolation method. For a true 
experimental spectrum, the calculation of r'(0) yields a set of 
compatible pairs of parameters, (c, £), and, if the concentra- 
tion is known independently, one determines £, as exemplified 
for the experimental data of Ref. [25j For an isotropic 3D en- 
semble of particles, the proposed criterium is shown to vanish, 
independently of the presence or absence of correlations. In 
this case, any significant deviation of r'(0) from zero would 



imply an anisotropy in the system. The knowledge of the sin- 
gle particle polarizability is required for the proposed anal- 
ysis, in particular, the particle radius is important. Here we 
considered metallic NPs described by the Drude model but 
the procedure can be applied also to semiconductor quantum 
dots with the polarizability determined by localized exciton 
resonances. 
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